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Abstract. In this paper we study subdiffusion in a system with a thin membrane. 
At the beginning, the random walk of a particle is considered in a system with a 
discrete time and space variable and then the probability describing the evolution of 
the particle’s position (Green’s function) is transformed into a continuous system. Two 
models are considered differing here from each other regarding the assumptions about 
how the particle is stopped or reflected by the membrane when the particle attempts 
to pass through the membrane fails. We show that for a system in which a membrane 
is partially permeable with respect to both its sides the Green’s functions obtained 
for both models within the continuous time random walk formalism are equivalent to 
each other and expressed by the functions presented in the paper: T. Kosztolowicz, 
Phys. Rev. E 91 , 022102 (2015), except the values defined at the membrane’s surfaces. 
These functions generate the boundary condition at the membrane which contains a 
term with the Riemann-Liouville fractional derivative which vanishes over sufficiently 
long time, this term is present even when normal diffusion occurs. This boundary 
condition can be interpreted as the particle’s passing through the thin membrane being 
a ‘long memory’ process for subdiffusion as well as for normal diffusion. We also show 
that for a system with a one-sidedly fully permeable membrane, the Green’s functions 
obtained within both models are not equivalent to each other and generate a boundary 
condition which does not contain the term with a fractional time derivative. 

We also present the generalized method of images which provides the Green’s 
functions for the membrane system, obtained in this paper. This method, which has a 
simple physical interpretation, is of a general nature and, in our opinion, can be used to 
find the Green’s functions for a system with a thin membrane in which various models 
of subdiffusion can be applied. As an example, we find the Green’s functions for the 
particular case of a ‘slow-subdiffusion’ process in a system with a thin membrane. 
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Normal diffusion or subdiffusion in a membrane system widely occurs in biology and 
engineering sciences (see, for example, HIEIE])- In this paper we consider subdiffusion 
in a one dimensional system with a thin membrane. 

The most commonly used definition of subdiffusion is that this process is the random 
walk of particles in which 


((Ax)2) = 




r(l + a)’ 

where ((Ax)^) is the mean square displacement of a particle, a is a subdiffusion 
parameter (subdiffusion exponent), is a subdiffusion coefficient, and T is the Gamma 
function. Subdiffusion occurs in media in which the particles’ movement is strongly 
hindered due to the complex internal structure of the medium, such as, for example, 
in porous media or gels HE]. The sub diffusion is often described by the following 
subdiffusion equation with the Riemann-Liouville fractional time derivative H 

dP{x,t;xo) _ d‘^P{x,t;xo) 

dt dx^ ' ^ ^ 

The Riemann-Liouville fractional derivative is defined as being valid for "d > 0 (here k 
is a natural number which fulfils k — 1 < -d < k) 

d^f{t) 1 


dp 


T{k - id) dP 


Jo 


( 3 ) 


There arise a problem how to set boundary conditions at a thin membrane (see, for 
example, H El H [ID]). 

The structure of Eq. ([3]) shows that the presence of the fractional derivative in 
a model can be interpreted as dealing with a ‘long-memory’ process. For a = 1, the 
process is usually identified as normal diffusion, which is considered as the Markov 
process. Normal diffusion can be interpreted as a particle’s random walk in which the 
mean square displacement of a single jump length and the mean frequency of jumps 
are both finite. However, there are processes in which the anomalously long waiting 
time for a particle to take its next step is entangled with the anomalously large length 
of jumps in a special way that provides a = 1 [6|. As concluded in [6] in order to 
define subdiffusion, relation ([1]) should be supplemented by an appropriate stochastic 
interpretation of the random walk process. Such a simple interpretation is given within 
the continuous time random walk (CTRW) model where the random walk is described 
by the probability density A(p) of a single jump length p and a probability density oj{t) 
of time r, which is needed for the particle to take its next step. It is assumed that for 
normal diffusion both distributions have hnite moments of all natural orders whereas 
for subdiffusion the mean value of uj{t) is infinite and the moments of \{p) are finite. In 
this paper we base our consideration on the random walk model on a lattice for which 
A(p) = |5(p — e) + \5{p + e) (at the vicinity of the membrane the definition of A(p) is 
slightly different), e is the distance between discrete sites; in this paper 5 denotes both 
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the Kronecker symbol (for a discrete space variable) or the delta-Dirac function (for a 
continuous space variable). The kind of diffusion process is dehned here by the function 
u;{t). Obtaining the results for a system with a discrete spatial variable, we pass into a 
continuous system in the limit e —?■ 0, using the formulas presented in this paper. 

In paper [10], subdiffusion in a system with a thin stopping membrane was 
considered; the membrane was assumed to be partially permeable with respect to 
both its sides, 0 < qi,q 2 < 1- The Green’s functions, obtained for this case, generate 
the boundary condition at the membrane which contains the term with the Riemann- 
Liouville fractional derivative; this term, which is present even when normal diffusion 
occurs, vanishes over sufficiently long time. This boundary condition can be interpreted 
as the particle’s passing through the thin membrane being a ‘long memory’ process. In 
this paper we extend the consideration to the case of a reflecting membrane, subdiffusion 
in a system with a one-sidedly fully permeable membrane is also studied within the 
models of reflecting or stopping membranes. 

There are three main aims of this paper. Firstly, we make a comparison between 
the Green’s functions obtained for the model in which the particle can be reflected from 
the membrane and the model in which the particle can be stopped by the membrane. 
We hnd the conditions under which these functions are not equivalent to each other. 
Secondly, we examine the boundary conditions at the thin membrane and check when 
these conditions do not contain ‘long memory’ terms. At this point we consider in detail 
the case of a one-sidedly fully permeable membrane. Thirdly, we hnd the generalized 
method of images which provides the Green’s functions for the membrane system derived 
in this paper within the GTRW model. Next, using the generalized method of images, 
we obtain the Green’s functions for a system with a thin membrane, in which the special 
case of ‘slow subdiffusion’ occurs. ‘Slow subdiffusion’ is understood here as a random 
walk for which all fractional moments of cu(r) are inhnite, (r^) = /q°° r^a;(r)(ir = oo for 

p>o nn. 

The organization of this paper is as follows. In Section |2l we present the general 
assumptions of the method. As an example, the Green’s function for a homogeneous 
system without a membrane is derived. In Section [3] we determine the Green’s functions 
for a system with a membrane which reflects particles with some probability. In Section 
m we show the Green’s functions for a system in which particles can be stopped by 
the membrane; the Green’s functions for a system with membrane which is partially 
permeable with respect to both its sides are taken from paper [lO]; the new result 
is obtained for the case of a one-sidedly fully permeable membrane. The generalized 
method of images is presented in Section El As an example of the usefulness of this 
method, we hnd the Green’s function for ‘slow-subdiffusion’ in a system with a thin 
membrane. A comparison of the Green’s functions obtained from the models presented 
in the Sections El and 0] and a discussion of the new results obtained in this paper are 
presented in Section El 
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We focus our attention on deriving the probability density (Green’s function) P{x, t] xo) 
of finding a particle at point x after time t under the condition that at the initial moment 
f = 0, the particle was at point xq (here we consider a particle’s random walk in a one¬ 
dimensional system). To hnd this function for a system with a thin membrane we use 
a particle’s random walk model in a system in which the time n and spatial variable m 
are discrete; after this, we pass from discrete to continuous variables. 

In a system with discrete variables, the particle’s random walk is described by a 
difference equation whose general form reads 

P„+i(m;mo) = '^Pm,m'Pn{rn'-.mo) , (4) 

m' 

where Pm,m' is the probability that a particle jumps from site m' directly to site 
m, Pn{m]mQ) is the probability of hnding a particle at site m after step n and mg 
is the initial position of the particle 


Poim^mo) = 5m,mo ■ 


(5) 


For normal diffusion and subdiffusion long jumps occur with a low probability, thus we 
assume that the particle can only jump to its adjacent position. It is not allowed to 
remain in the current occupied position in the particle after the time at which the jump 
should be executed unless the particle is stopped by the membrane. Such a model is 
useful in modeling normal or anomalous diffusion processes in a system with a partially 
permeable thin membrane. The reason is that Eq. (jl]), applied to describe a random 
walk in a membrane system in which the homogeneity of the system is impaired at a 
single point, is solvable by means of the generating function method |l2l [131 [H] . After 
solving these equations (more precisely, after determining the generating function for 
the equations) one can make the transition from discrete to continuous variables. The 
transition from discrete to continuous time is performed as in the CTRW method while 
passing from a discrete to continuous space variable is performed using the formulas 
presented in this paper. The reason for the introduction of such formulas is the specihc 
behavior of the particles in the vicinity of the partially permeable membrane. 

In [15] the random walk in a homogeneous system was considered assuming 


- Ia 

Pm,m' 2^5m—l,m' 4 “ 2^5m+l,'i 


( 6 ) 


The considerations were also extended to the case of normal diffusion in a system with a 
fully reflecting or fully absorbing membrane. In these systems, the ‘parity rule’, meaning 
that the number of steps n and the distance traveled by a particle |m—mo|, are both even 
numbers or both odd numbers, is fulhlled. For the system with a fully reflecting or fully 
absorbing membrane the Green’s functions were obtained using symmetry arguments. It 
was shown [T5| that if the wall were removed, for each particle’s trajectory which passes 
the line representing the position of the wall, there exists a symmetrical trajectory 
with respect to the wall, which is solely located in the half-space bounded by the wall. 
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in which the particle started its random walk. In this way it was shown that the 
Green’s functions for these systems is the sum (for the system with a fully reflecting 
wall) or difference (for the system with a fully absorbing wall) of the Green’s functions 
obtained for the homogeneous system without the wall, and the initial points Xq of these 
functions are located symmetrically with respect to the wall. The Green’s function can 
be interpreted as a normalized concentration profile of particles initially located at the 
initial point, which means the concentration of N particles, divided by iV, iV 3> 1. 
By differently describing the method of images, the idea of this method is that the 
wall is replaced by the additional particles’ source located in such a way that the 
particles’ concentration generated by the particles’ source located at mo and by the 
additional particles’ source both fulfil the appropriate boundary condition at the wall: 
the vanishing of the particles’ flux at the reflecting wall and the vanishing of the particles’ 
concentration at the absorbing wall. The method of images appears to be useful tool in 
deriving Green’s function for a system containing two walls. 

The situation is different when a thin partially permeable membrane is located in 
the system. The presence of the membrane can be taken into account in the model in 
two different ways which depend on assumptions about how a particle is stopped by the 
membrane. Assuming that a thin membrane is located between the N and the + 1 
site, there are two possibilities concerning the influence of the membrane on a particle’s 
movement. In the first model, the particle can be ‘reflected’ by the membrane, which 
means that after the jump the particle always changes position. If the particle’s attempt 
to pass through the membrane from the N to the + 1 site fails with probability gi, 
the particle is reflected by the membrane from the N to the A^ — 1 site; if the particle 
tries to jump through the membrane from the A^ + 1 to the N site and this attempt 
fails with probability q 2 , then the particle is reflected from the A^ + 1 to the A^ + 2 site. 
In this case we obtain (see Fig. 1). 



Figure 1. System with reflecting membrane, more detailed description is in the text. 
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Pm,N+l 


Pm,m' 


1 + gi 
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1 — Q'2 r 1 + Q'2 r 

- - -d- - - 0 

1 1 
l,m' T 


V+2,m ) 

, m' N, N + 1 . 


For this model the ‘parity rule’ is satisfied. 


( 7 ) 

( 8 ) 
( 9 ) 
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For the model in which the particle can be stopped by the membrane, the particle 
remains in its position after the ‘jump’ when the passing through the membrane has 
failed with probabilities qi and q 2 . In this case we have (see Fig. 2). 



Figure 2. System with stopping membrane, more detailed description is in the text. 


Pm,N 

Pm,N+l 

Pm,m' 


^ A I A I ^ A 

-ON-l,m + —ON,m H- - -OV+l.m , 

-^0N+2,m + -^Ojy+i^rn H-^- Ojy^rn , 


( 10 ) 

( 11 ) 

( 12 ) 


For this case the ‘parity rule’ is not satished. In the following, we will check if the 
assumptions involved in the models provide any noticeable differences between the 
Green’s functions for a continuous system. 

To solve the differential equations we use the generating function which is dehned 


as 


S(m, z; mo) = y) z''P„(m; mo) ■ 


(13) 


n=0 


To pass from a discrete to continuous time we use the following formula 

00 

F(m,f;mo) = J2Pn{m;mo)^n{t) , (14) 

n=0 

where <Fn(^) is the probability that the particle takes n jumps over time interval (0,f). 
This function fulhls the following relation m ^n{t) = !lu{t - t')Qn{t')dt\ where 
U{t) = 1 — jQu{t')dd is the probability that the particle does not perform any jump 
over time interval (0,f), u}{t) is the probability distribution of the time which is needed 
for the particle to take its next step. Quit) denotes the probability that the particle 
performs n steps over this time interval and the last step is taken exactly at time f. 
Quit) = — t')u{t')dt', n > 1, Qi{t) = Qo{t) = 6{t). In terms of the 

Laplace transform, C[f(t)] = f{s) = f(t)dt, the function ^n(t) reads [ 1 ] 

--[^(s)] • (15) 


$ 
^ r 


s = 


Combining the Laplace transform of Eq. flTTll with Eqs. (|T^ and (|T5il we obtain 


P(m, s; mo) = 


1 — a)(s) 


S (m, d;(s); mo) 


(16) 








Subdiffusive random walk in a membrane system. 

To pass from a discrete to continuous spatial variable we suppose 


7 


X = em , Xo = emo , xjsi = eN , (17) 

where e denotes the distance between discrete sites. To pass form probability P{m, s; mo) 
to spatial probability denisty P{x, s; Xq) we use the following relation valid for small 
values of e 


F(m, t; mo) = eP(x, f; Xq) , (18) 

and hnally we take limit e —> 0. Further considerations are performed assuming that 
s is small, which corresponds to the case of large time due to Tauberian theorems. In 
practice, the limit of small s means that only the leading terms with respect to this 
variable will be present in the Laplace transform of the Green’s function whereas the 
limit of small e means that this parameter will be absent in this function. 

As an example we consider subdiffusion in a homogeneous system without bias. 
Starting from the following differential equation 


F„+i(m;mo) = ^F„(m - l;mo) + ^F„(m + l;mo), 
from Eqs. ([5]), flTdll and flT^ we get (see [I2l [16] and the Appendix in this paper) 


S(m, z; mo) = 


\/l — 


where 


1(0 = 


vr 


(19) 


( 20 ) 


( 21 ) 


In the standard continuous time random walk model it is assumed that a;(s) for small s 
reads |1] 


a;(s) = 1 — Tr,s^ 


( 22 ) 


where is a parameter which, together with a, fully characterizes time distribution 
u{t). Eqs. fl2Tl) and fl22l) provide, for small values of s 

T]{uj{s)) = 1 - \/2raS“ . (23) 


The subdiffusion coefficient is dehned to be 


= 


2r„ ■ 

From Eqs. (IT6l) - fl2T|) we observe 


F(x, s;xo) = 




5a/2 X 


U-^nl 


(24) 


(25) 


Formally, the transition to a continuous variable was obtained by calculating limit 
e —> 0. However, due to Eq. fl2TD this provides Tq —> 0 which means that (xit) 
is beyond any physical interpretation. In order to avoid problems of interpretation, we 
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assume that e is finite but small enough so that the following relation is satisfied (see 
also the discussion presented in [TO] ) 




1 -e- 


Y 


n-foi 




a/2 


\/Da J 

From Eqs. fl2S]) and 0261) we have 

^ 1 
P{x, s;xo) = 


-0 V Da 


(26) 


(27) 


The inverse Laplace transform of Eq. (|26|) . calculated by means of the following formula 

im 


c 


-1 




1 oo ^ 

= ^ g kW(-k/3 - 0 V ’ 


(28) 


a, /5 > 0, reads 

P{x, t; Xo) = /a/2-1,a/2 . (29) 

We add that the function can be also expressed in terms of the Fox function. 
Function fl29|l is the solution of a fractional subdiffusion equation with a Riemann- 
Liouville fractional derivative ([2]) for the initial condition P(a:,0;a:o) = d{x — xq). 


3. Random walk model with particle reflection from a membrane 

Let us suppose that the thin membrane is located between the N and the iV + 1 site. 
We assume that the particle can be reflected from the membrane with probability qi 
when trying to jump from the N to the iV + 1 site and with probability q 2 when trying 
to jump from the iV + 1 to the N site (see Fig. [1]). Combining Eq. (]1]) and Eqs. ( 0-0 


we can observe 

Pn+i{N - l;mo) = ^PniN - 2; mo) + ^—^Pn{N;mo), (30) 

F„+i(A^;mo) = ^Pn{N - 1; mo) + Pn{N + 1; mo), (31) 

F„+i(A^ + l;mo) = ^F„(iV;mo) + ^Pn{N + 2; mo), (32) 

Pn+iiN + 2; mo) = — ^^Pn{N + 1; mo) + —Pn{N + 3; mo), (33) 

F„+i(m;mo) = ^F„(m - 1;mo) + ^F„(m + 1;mo), (34) 

m 7 ^ — 1, A^, + 1, + 2, 

the initial condition is given by Eq. ([5]). 


In the following, the functions S and F will be labelled by the indexes ij, which 
denote the signs of m — A^ and mo — N, respectively. We consider the case of mo < N. 
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From Eqs. (E]), ffT^ and (15U]) - 
in the Appendix) 


we obtain (the details of the calcnlation are presented 


T 

1 

1 

73 


S — {N,z]mo) 


S+-{m,z]mo) 

^m-mo 

S+-{N + 1 , 2 ;; mo) 

„V+1—mo 


2N—m—rao 


\/l — 


, m < N — 


where 


_ di- q2rf{z) _ 1 - q2Tf{z) 

d^i\z) — - - irrz ) Aijv(^) — 


K2{z) = 


1 - qiq2yf{z) 
(l-gi)(l + g2) 


, K2n{z) — 


1 - qiq2rf{,z) ’ 
1 - gi 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


1 - qiq2r]‘^{z) 1 - qiq2'n‘^{z) 

Functions fl3^ and fl40|) play a key role in the Green’s function for the membrane 
system, because only these functions depend on the permeability coefficients of the 
membrane. The main difficulty of finding the Green’s functions for the continuous 
system is to hnd a suitable form of functions fl39l) and fHOl) for e —)■ 0. We assume 
that in this limit, the Green’s function must depend on the parameters of membrane 
permeability. As we will see below, the case of gi, g 2 > 0 and the case of gi = 0 or g 2 = 0 
should be considered separately. 

Similarly to [10], the probabilities characterizing membrane permeability are 
assumed to be functions of e; this function for small values of argument e reads 


gi(e) = 1 

52(e) = 1 


71 
e'" 

72 


(41) 

(42) 


7 i and 72 being the membrane permeability coefficients defined for the continuous 
system, and a being the parameter to be determined. The reason for the introduction 
of these equations is that the frequency of jumps performed by the particle increases 
to inhnity when e drops to zero. For very small value of e, the frequency takes an 
‘anomalously’ large value |T 0 ]. A very large number of attempts to pass through the 
partially permeable membrane made over an arbitrarily small time interval means that 
the particle passes trough the membrane with probability equals to one; then, the 
membrane loses its selective property. 
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K,{u{s)) = 


Kin{uj{s)) = 


K2{u:{s)) = 


^a-l 1 

^---1 
71 12 J 

1 + 2^ 

j sa/2 
/ 

^a-l 1 

f+ + +') 

171 12) 

1 + 2^ 

Ua.l1 
' Da 


^CT—1 1 I 9 / <5^^^ 


t-cr— 1 


(^ + i) + 2 ^ 57 - 


2g<7-li 

71 


■‘(7 + i) + 2v^’ 




^(T— 1 J_ 
71 


(43) 

(44) 

(46) 

(46) 


(^+ +2v'^ 

The only possibility that Green’s fnnctions depend on the parameters of membrane 
permeability in the limit of small e, also for the case of symmetric membrane for which 
7i = 72 , is cr = 1. Then, form Eqs. flT3|) -fH6 |l we get for small s 

4., „a/2 

(47) 

(48) 


m \\ 71-72 , 

i7i(ci;(s)) =-^-h 


71 + 72 (7i + 72)\/A^’ 


rm \\ 7i , 2y2y^s"+ 


Muis)) = 


7i+72 (7i + 72)\/A^’ 

2y2 4y2y^s"+ 


K2n{uj{s)) = 


7i + 72 (7i + 72)\/A^’ 

72 2y2y^s“+ 


71+72 (7i+72)\/A^’ 


(49) 

(50) 


wli6r6 'y = 
wiieie 7^ 71 + 72 * 


From Eqs. flTB]) - f[TSD . fl5^ - (l55]l and flT7|) - fl5(I]) we can observe 


P—{x,s;xo) = 


27d;^si-“/2 


lx —Xq I 

0 V Dql 


( 61 ) 


i^7i — y2 4y2y«,s“+ \ (2xj^-x-x^)s°‘/‘^ 

-1- 1= ) e 'JDa 

v7i + 72 (7i+72)v7^/ 


, X < Xn 


b ( \ ^ 2y2y^s"+ . 

P—(xjv,g;a(o) = ^ 1 , ., + -—^^^^(15?) 


27d::s1-“/2 


^71 + 72 (7i + 72)\/A^4 


P+_(x,s;xo) = 
X > Xn 1 


_ 272 _ 4y2y^s°+ 

2\/l7)s^-"/2 _y^ + y2 (y^ + '32)\fD^_ 


P+-{xn,s-,Xo) = 


2^/P;sl-“/2 


72 2y2y^s“+ 


71+72 (7i+72)\/^ 


e v73^ , (53) 


{x]^-Xf^)s°‘/‘^ 

e vda (54() 
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The inverse Laplace transform of the above functions read 
P - {x,t]Xo) = /a/2-1,0/2 \t',^ 


7i - 72 


2v^(7i + 72) 
2727 . 


2^ 

/a/2-1,a/2 t] 


\/ Da } 

2xn — X — Xq' 


-Da (71 + 72 ) 
P—{xN,t]Xo) = 


w r . 2 xn - x-xo\ 

/a—l,a/2 I I , X Xj^ , 


\/Da 


7i 


2v^(7i + 72) 

Xn — Xq 


/a/2-1,a/2 t] 


Xn - Xq 




P+-{x,t]Xo) = 
2'l2lw 


72 


\/D^(7i + 72 

^ X — Xq 


/a/2-1,a/2 t] 


Da (71 + 72 ) 
P+-{xN,t-,XQ) = 
l2lw 


fa- 1 ,a /2 t] 


\/D^ / 




3^ - 3^0 \ 

x/D; ) 

X > Xn , 


72 


2VD^(7i + 72) 

Xn — Xq 


/a/2-1,a/2 t] 


Xn - Xq 


fa- 1 ,a /2 t] 


( 66 ) 


( 66 ) 


(57) 


(58) 


£'»(7i + 72)'’“^“''7'’ ) ' 

It is easy to see that the Green’s functions lose their continuity at membrane surfaces 

1 


P±_{xN,t;xQ) = 2^±-ix 


■xff,t;xQ). (59) 

When the membrane is one-sidedly fully permeable, e.g. in the case of gi = 0 or 
q 2 = 0, the dependence of functions fl3^ and fITOl) on the parameter e is different than 
in the previous case, which we will show in the next sections. 


3.2. The case qi = 0 and 0 < ^2 < 1 

Supposing gi = 0 and g 2 > 0, combining Eqs. 
obtain over a limit of small values of e and s 


, and 


, and 


Ki{u{s)) — -1 H-h 2eJj— , KiN{dj{s)) — -h 2ew— 

72 V Da 72 V Do 


we 

(60) 

(61) 


A’ 2 (w(s)) ^2-, A' 2 iv(w(s)) = 1 . 

72 

The Green’s functions are dependent on parameter 72 over the limit of small values of 
e only if a = 0; in this case, the functions are independent of parameter s. Proceeding 
as in Sec. 13.11 we obtain 

^ |a: —a:o|\ 


P - {X,t-,XQ) — fa/ 2 -l,a /2 t] 


^ ~ 7 V ’ 


V Da j 

2Xn — X — Xq' 

\/^ 


( 62 ) 


, X < Xtv , 
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P—{xN,t-, Xo) 


1 . ( Xn - Xo 


P+-{x,t;xo) 



/o/2-l,a/2 



3^ 

^/^ ) 


X > Xn , 


P+_{xN,t;xo) 


( xn- Xo\ 


(63) 

(64) 

(65) 


3.3. The case 0 < gi < 1 and q 2 = 0 

Supposing q 2 = 0 and qi > 0, from Eqs. (|2^ . (j2^ . ([39]), (HUjl . and (HTjl we obtain 

Ki{u{s)) = 1 - — , iriAr(w(s)) = 1 , (66) 

7i 

K2{u:{s)) = K2N{do{s)) = - . (67) 

7i 

Thus, Eqs. (I66|) and (167|) depend on the membrane permeability coefficient over the 
limit of small values of e only if cr = 0. After calculations, we have the following Green’s 
functions 


P - {x,t;Xo) — /a/2-1,a/2 ( t; 


\x — Xo\ 


, /1 1 W (^ 2xn - X - xo\ ^ ^ ^ 

[ tJ r 7m )' ' 


P - {xN,t',Xo) — /a/2-1,a/2 ( t; 


Xn - Xo 




n / N 1 ^ ( X — Xo 

H— (3') t, Xo) Jal2—l,al2 I t, . I , X Xn ; 

7l V v Pa J 




P+-{XN,t\Xc,) = Pa/2-1,a/2 (*’ 


( 68 ) 


(69) 

(70) 

(71) 


4. Random walk model with particle stopped by a membrane 

In |T0] we considered a particle’s random walk model in a system with a thin membrane, 
in which the particle trying to pass the membrane from the N to the 77 + 1 site, may 
be stopped by the membrane with probability qi or pass through the membrane with 
probability 1 — gi. The particle’s halting by the membrane means that the particle 
does not change its position after its ‘jump’. When the molecule is trying to jump 
from position + 1 to 77, the probability of the blocking of the particle through the 
membrane is q 2 , and the probability of passage through the membrane equals 1 — g 2 - 
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From Eqs. (|1]) and (fTU]) - fll 2 p we have the following difference equations 

P„+i(m;mo) = ^Pn{m - l;mo) + ^Pn{m + l;mo), N,N + |72) 

Pn+i{N; mo) = ^Pn{N - 1; "^o) + ^—^Pn{N + 1; mo) 


+ |F„(iV;mo), 

Pn+iiN + 1; mo) = —mo) + -Pn{N + 2; mo) 
+ —Pn{dd + l;mo). 

The generating functions for the equations fl72p - fl74p read ra 


S _(m, z; mo) = 




VT^ 


+ ^ 1 ( 2 ^) 




S'+_(m, 2 ;;mo) = A 2 {z) 




m—mo 


» 




where 


Ai(z) = 
A 2 {z) = 


qi - q2v{z) 


1 - (gi + g 2 - ’ 

{1 + ri{z)){l - qi) 


1 - (gi + g2 - i)h(^) ’ 
4 . 1 . The case 0 < gi < 1 and 0 < ^2 < 1 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 


In the following 71 and 72 denote the permeability coefficients of the stopping membrane, 
dehned for a system with a continuous spatial variable. As shown in [TU], the 
probabilities of passing through the membrane should be chosen as the following 
functions of e 


gi(e) = 1--^ , q 2 {e) = 1 - ^ , 
7i 72 


which, together with Eqs. 


, (1771) . flTSp . and (17^ provide 


\ \\ 71 -72 , 2727^ s “/2 

Ai(a;(s)) = ' . , + 


71 + 72 (7i+72)\/7^’ 

272 _ ‘^l2lwS°‘/‘^ 

71+72 (7i+72)\/^’ 

where 7 ^ = . The Green’s functions read m 


A 2 (u(s)) = 


P__(x,t;xo) = 


2 V^ 


/a/ 2 - 1 ,a/2 t] 


|X — Xol 

V Da 


(79) 

(80) 
(81) 

( 82 ) 


, 7i - 72 . f, 2 xn -x-Xo 

A- ---^/a/2-l,a/2 t] 

7l + 72 V 

, l 2 lw 


(71 + 72)-D, 


■/a—l,a/2 I A 


\/Da j 
2xjq — X — Xq 
VP a 
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P+-{x,t;xo) = 


272 


7i + 72 2 ^/ D ^ 


1 r ( X-Xq' 

Jal2-l,al2 1 


\/Da ) 


l2ln 


(71 +72)^0 


/o—l,a/2 I t, 


x-xo \ 

) 


(83) 


The Green’s functions Eqs. (18^ and (l83|) coincide with the Green’s functions 
obtained for the system with a reflecting membrane, Eqs. fl55ll and flHTll . respectively, if 
7 i = 71/2 and 72 = 72 / 2 . 

As in Section [3l the case of a membrane which is one-sidedly fully permeable should 
be considered separately. 


4 . 2 . The case qi = 0 and 0 < q 2 < I 
We suppose that qi = 0 and 


q2{e) = 1 - — 
72 


Form Eqs. 


, dZID, (Hi and 


we obtain 


Ai(z) = 


72 ' 

1 + 

72 


, A2(^) = 


2 -e, 


1 + 

72 


(84) 


(85) 


As previously, functions Eq. flHj depend on the parameter of membrane permeability 
only for a = 0. In this case, the Green’s functions read 


P—{x,t;xo) = 


1 


2V^ 


fal 2 -l,al 2 t] 


\X — Xq\ 

y/Da 


72-1 f 2 xn - X -xo' 

jQ/2-l,a/2 t, ■ 


72 + 1 


P+-{x,t-,xo) = 


'/Da 


72 


(72 + ^)/T/ 


/o/2-l,a/2 


x--^\ 

/^a ) 


( 86 ) 


(87) 


4-3. The case 0 < gi < 1 and q 2 = 0 

Supposing g 2 = 0 and 

gi(e) = 1 - ^ , 

7i 

proceeding as in Sec. 14.21 we obtain 


1 - e‘"A 2e^A 

Ai(A(5)) = , A2(+(.)) = 


1 + e'’ 


71 


1 + e'^A 

71 


( 88 ) 


(89) 


Similarly to the previous case, we have to take a = 0. The Green’s functions read 


P—{x,t;xQ) = 


1 


2 /^ 


/a/ 2 - 1 ,a/2 t; 


\X — Xq\ 
/Da 


7i - 1 2xn - x-Xo 

+ - -^Ja/2-l,a/2 ' 


7l + 1‘ 


V Da 


( 90 ) 
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P+-{x,t;xo) 


1 

(71 + 1 )\/-Dq 


L 


12 - 1 ,al2 



X - Xo \ 

vd: } 


(91) 


5. Generalized method of images 


As we briefly discussed in the Introduction, the classical version of the method of 
images has been used to determine the Green’s function in a system containing a fully 
reflecting or fully absorbing wall. In this method, Green’s function is considered to 
be the normalized concentration of diffusing particles beginning their movement from 
initial point xq at t = 0 . The main idea of this method is to replace the wall by an 
additional source of particles in such a way that the concentration of particles generated 
by all particle sources occurring in the system, fulfils the boundary condition which is 
assumed at the wall. For a system with a fully reflecting wall, the boundary condition 
at the wall is that the diffusive flux vanishing at the wall, J{xN,t]Xo) = 0 , xn is the 
position of the wall. The same effect would be achieved if the wall were to be replaced 
by an additional source of particles located symmetrically to point Xq with respect 
to the wall. For a system with a fully absorbing wall, the boundary condition reads 
P{xp^,t-, Xq) = 0 , the additional source located symmetrically to point xq with respect 
to the wall should be substracted from the Green’s function representing the particles’ 
source located at Xq. Thus, assuming that Xq < xat, the Green’s function for the above 
mentioned cases can be written in the following compact form 

P{x, t; xo) = Po(3^, t; Xq) + <;Po{x, t] 2xn - Xq) , (92) 

where ^ = 1 is for a fully reflecting wall and <^ = — 1 is for a fully absorbing wall, Pq 

denotes here the Green’s function for homogeneous system with removed wall. 

The form of the Green’s functions, presented in Secs. |3] and IU shows that these 
functions can also be determined using the method of images which is understood here 
as a replacement of the membrane by the additional source function. Analyzing the 
structure of the Green’s functions obtained in the previous section we note that the 
functions can be expressed by the following equations 

P—(x,t;xo) = Po(x,t;xo) + Pc(x,t;2x]v - Xq) , (93) 

P+-(x,t;xo) = Po(x,t;xo) - Pc(x,t;xo) , (94) 


where Pc denotes a ‘compound’ source function. This function has the following 
structure 


Pc(x,t-,xa) = Ki,P(,(x,t;X(,) + KcPc(x,t',Xa) , 


(95) 


where 


Pcix.t-.XQ) 


—Po{x,t-,Xo) 


(96) 


is a ‘gradient source function’. Parameters Kq and kq depend on the membrane 
permeability coefficients. In general, kq can be interpreted as the relative measure of 
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the asymmetry of the membrane; /to = 0 for the symmetrical membrane. The particular 
cases are the following: for 0 < qi,q 2 < 1 for the system with a reflecting membrane we 
observe 

7i - 72 47 i 

Ko = -^ 

71+72 (1 + 71/72) 

For the system with a stopping membrane, Eqs. (|Mll - (177)l are still valid, with 
7 i ,2 = 71,2/2. For the case of a fully reflecting wall we have 71 — > 00 , then /to = 1 

and Kg = 0 which means that a function P _takes the form of function fl9^ with 

= 1 and P+_ equal to zero. For a partially absorbing wall there is 72 —> 00 and 
0 < 7 i < 00 which gives Kq = —1 and kq = 471 . When the membrane is one-sidedly 
fully permeable, then kq = 0 for all cases, but parameter kq is not the same for both the 
models mentioned above. For the model of the reflecting membrane we have: for qi = 0, 
g 2 > 0 there is kq = —I + I /72 and for ^2 = 0, gi > 0 there is kq = 1 — 1 / 71 . For the model 
of the stopping membrane we have: for gi = 0, ^2 > 0 there is Kq = (—l+l/ 72 )/(l+l/ 72 ) 
and for ^2 = 0, gi > 0 there is /to = (1 — l/ 7 i)/(l + l/ 7 i)- 

We assume that the method of images may be used to determine the Green’s 
function for various models of subdiffusion in a membrane system. As an example, we 
consider the ‘slow subdiffusion’ model, in which the waiting time for the jump is given 


by the following function u{t) = 1 — 
f S> /i, reads 


In/i 




r—1 


Inf J 

Let the subdiffusion coefficient be defined 


, p, r > 1 [ID, which, for a long time. 


(98) 


the Green’s function for these processes for a homogeneous system without a membrane 
reads [11] 

|a;-xo| 

I T 

Fo((c,f;xo) = 


\/^ 


aint)’’ 


2GF>r,/.(hif) 


( 100 ) 


This Green’s function provides the relation ((Ax)^) = 2Dr^^{\atY ^ [TT|. Using Eq. 
we have 




PG{x,t]Xo) = 




( 101 ) 


The Green’s functions for the slow subdiffusion process in a system with a thin 
membrane, derived by means of the generalized method of images, read 


F__(a;,f;xo) = 


Kq 


2 i/F, 


Kq 


(lnf)'’“i 

-^ 0 ) 


r,/! 
(2xjV 


\/hrU 


(Int)’- 


2uFr,/i(lnf) 


r—1 




+ 






{2xj^—x 


I^G 


^ 0 ) 


2Dr,Y'^nt) 


r—1 


\/^ 


aint)’’ 


( 102 ) 
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P+-{x,t;xo) = 


I- Ko 




Kg 


2Dr,i,{\nty-^ 


2 sJ Dr,y\nty-^ 




(103) 


coefficients kq and kq are defined by Eq. (1771) for a both-sided partially permeable 
membrane, and are defined by the equations presented in Sec. [5]just below Eq. (177)1 for 
a one-sidedly fully permeable membrane. 

The following question arises: whether the generalized method of images provides 
functions for the slow subdiffusion that is consistent with the functions obtained by 
means of the random walk model presented in this paper? To have the answer to this 
question, which will appear to be positive, we perform the following consideration. From 
Eqs. (El]), (|M1) and the Laplace transform of Eq. (IM|) . which, due to Tauberian theorem, 
for small values of s reads P 


a;(s) 


1 - 


( \ 
l,ln(l/s); 


(104) 


we obtain over the limit of small values of e and s 
r]{u{s)) = l-e ^ 


Dr,yin{l/s)y 1 


(105) 


From Eqs. (l3^ - (IT2l) . (TTTll - dTT) . (11041) and (llOSp . after simple calculation, we obtain 
functions iLi,iv, 2 , 27 v (for the model of the reflecting membrane) and Ai ^2 (for the 
model of the stopping membrane) for the ‘slow subdiffusion’ case, which takes the 
form of the analogical functions presented in Sec. [3] and Sec. 0] after substituting 

Further calculation, performed using Eqs. (fT6|) . (fT7|) . (1T8|) . 


,yDrAin{i/s)r 

(I104p . (I105p . and functions i^i,iAr, 2 , 2 Af and A 12 described above, provides the 
Laplace transforms of Eqs. (I102p and (I103p . 

Figures |3] and m are merely on illustration of the Green’s functions but show some 
general regularities of subdiffusion in a membrane system. Membrane permeability 
property is less manifest in a system in which slow subdiffusion occurs. This is due to 
the fact that when a particle tries to jump through a thin membrane, the probability 
of a particle being stopped by the membrane is the same as for the different processes 
occurring in the membrane system but a subsequent attempt at a jump is expected to 
be much later in the system in which slow subdiffusion occurs. So, the probability of 
a particle’s passing through the membrane over a certain time interval is less for this 
process compared to the ‘standard’ subdiffusion process. The qualitative comparison 
of the changes of the Green’s functions for different orders of the time variable is also 
interested. For ‘slow subdiffusion’ the changes are relatively small (see Fig. 0]), whereas 
for ‘standard’ subdiffusion the scale of the changes is much higher (see Fig. [3]). 
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X 


Figure 3. Plots of the Green’s functions for subdiffusion occurring in the system with 
a thin membrane (1551) and ((571) . obtained for a = 0.9, = 0 . 001 , 71 = 0 . 8 , 72 = 0.3, 

xq = —0.5 for times given in the legend, the thin membrane, located at xjv = 0, is 
represented by the dotted vertical line, all quantities are given in arbitrary chosen 
units. 



X 

Figure 4. Plots of the Green’s functions for ‘slow subdiffusion’ occurring in the 
system with a thin membrane (|102l) and (jl03|l . here r = 1.9, Hr./x = 0.001, the other 
parameters are the same as in Fig. [3l 
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In this study we consider subdiffusion in two systems with a thin membrane. The 
differences between the systems are related to the assumption concerning the process of 
the reflecting or stopping of a particle by the membrane when an attempt to pass the 
particle through the membrane fails. In the first model, the particle can be reflected off 
by the membrane with a certain probability, the second model assumes that the particle 
can be stopped, with a certain probability, at the membrane. These assumptions can be 
associated with the following physical interpretation. If there is a slight repulsive force, 
generated by the membrane and acting on the particle, then the model of the reflecting 
membrane can be used. Otherwise, including the case in which a small attraction of 
the particle exerted by the membrane is assumed, the stopping membrane model can 
be used to describe subdiffusion in a membrane system. 

The differences between the models are manifested in two cases: 

(i) The probability of finding a particle at the surface of the membrane is different 
for both models when the membrane is partially permeable with respect to both 
its sides, i.e. when qi,q 2 > 0. The Green’s functions are discontinuous at the 
membrane surface for the model with a reflecting membrane (see Eq. (15^ ) and 
continuous for the model with a stopping membrane. 

(ii) For the system in which the membrane is fully one-sidedly permeable, the Green’s 
functions for both models take a form which can be expressed by the function 
obtained within the generalized method of images fl93|) - (l95D if the ‘gradient source 
function’ Pq is omitted. However, the Green’s functions obtained for the model 
with a reflecting membrane and the ones obtained for the model with a stopping 
membrane are not equivalent to each other. The reason being that if is not possible 
to find the relation between the parameters 72 and 72 (for qi = 0 ) in such a way 
that Eqs. fl62|) and (l 6 T|) coincide with Eqs. fl 86 D and (IHTD . respectively. A similar 
remark can be made for the case q 2 = 0 . 

For the case qi,q 2 > 0, the Green’s function for both models, Eqs. (l55il and (1571) 
for the model with a reflecting membrane and Eqs. fl82|) and (l 8 ^ for the model with a 
stopping membrane, coincide with each other, respectively, in region (— 00 , XAr)U(xAr, cxd) 
if 7 i = 7 i /2 and 72 = 72 / 2 . Since the membrane permeability coefficients are expected 
to be determined from experimental data when the model is used to describe subdiffusion 
in a real membrane system, it does not matter which of the models will be used in the 
modelling unless the probability of finding diffusing particles at the membrane surfaces 
is not considered. 

To study the subdiffusion in a membrane system we use the model with a 
discrete time and space variable, next we transform the Green’s functions to continuous 
variables by means of the formulas presented in this paper. Such a model seems to be 
oversimplified. However, for the homogeneous system without a membrane, it provides 
results which can be derived by means of more ‘realistic’ models of subdiffusion. Thus, 
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we assume that the model used in this paper provides the Green’s functions which are 
useful in the modelling of subdiffusion in a system with a thin membrane. 

The method of images appears to be useful tool for determining the Green’s 
functions in a membrane system for a various model of subdiffusion for which the waiting 
time probability density is given by Eq. (EH]). It was shown that the generalized method 
of images allows us to obtain the Green’s function for the membrane system for a ‘slow 
subdiffusion’ process. We hypothesize that the generalized method of images, presented 
in this paper, is applicable to other subdiffusion models. 

The considerations presented in this paper concern the case of Xq < xp^. In 
general, the coefficient 71 can be dehned as a coefficient which controls the membrane 
permeability when a particle tries to pass through the membrane from a region of the 
initial particle’s location to the opposite region, and 72 is a membrane permeability 
coefficient when the particle moves in the opposite direction. In [10] it was shown that 
for the case of gi, ^2 > 0 the Green’s function for the model with a stopping membrane, 
Eqs. fl82l) and dESj), fulhl the following boundary condition at the membrane 

P—{xkp,t;Xo) = XiP+_{xjp,t;Xo) + X2- -’ ( 106 ) 


where Ai = 5^ and A 2 = -4k=- This boundary condition is complemented by the 
condition of flux continuity at the membrane. The considerations presented in this 
paper show that boundary condition Eq. fllObh is fulhlled by the Green’s functions at 
the reflecting membrane for gi, ^2 > 0 if we assume 71 = 271 and 72 = 272 . The second 
term in the right-hand side of Eq. (Iin 6 |l , which vanishes over a sufficiently long time (see 
the discussion presented in [lO]), represents the ‘additional’ memory effect created by 
the membrane (despite the fact that subdiffusion is a long-memory process itself). As 
mentioned earlier, for the case of gi = 0 or g2 = 0, the Green’s functions take the form of 
Eqs. (I93|)-(l95|) with kq = 0 for both models. It is easy to see that these functions fulhl 
the boundary condition (llObh in which A 2 = 0. Thus, the ‘additional memory effect’ is 
not created by the one-sidedly fully permeable membrane. 

For the case of xq > xn, the functions and boundary conditions can be obtained 
from the functions presented in this paper when, due to symmetry arguments, the 
following conversion is made: 71 —> 72 , 72 —> 7 i (or 71 —> 72 , 72 —> 7 i), 
xo —> xo — X, X — xtv — > xn — X and xn — xq —> xq — xn- 


X 
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Appendix 

Using Eqs. ([5]), ffT^ . and we obtain 

z z 

S{m,z;mo) - 6m,mo =-S'(m - 1, z; mo) +-S'(m + 1, z; mo), (107) 

m N — 1, N, N + 1, N + 2, 


S{N - l,z-, mo) - 6N-i,mo = ^S{N -2, z] mo) + A, z; m()ip8) 


S{N,z;mo) - 5 n, 


mo 


= -l,z-,mo) 


(109) 


- q2) 


S{N + l,z;mo), 


S{N + l,z-,mo) - dN+i,mo = ^^^-y-^5'(A,z;mo) + |S'(A + 2, z; m()i).0) 


S{N + 2, z] mo) - SN+2,mo = A + 1, mo) 

+ ^S{N + 3, 2 ;; mo). 


(Ill) 


To solve Eqs. (I107p - (llll|l we nse the following generating fnnction with respect to 
the space variable 


G(-u, mo) = ^ u^S(m, z]mo) . 


( 112 ) 


The generating fnnction S can be obtained by means of the following formnla 

1 r G{u, z] mo) 


S'(m, 2 ;;mo) = ^ / 
27rt Jk 


,m+l 


-du , 


(113) 


IK{0,1) U'- 

where integration is carried ont along the nnit circle K centered at point 0 = (0, 0) in 
order to be consistent with an increasing argnment of a complex nnmber. From Eqs. 
fll07p - flll2p we obtain 

,mo {u^~^ — M^+^) 


G'(M,z;mo) = 


u 


1 - - 
2 


-S{N + 1 , 2 ;; mo) 

Using the integral formnla 
1 


(“+h 

zq2 {u^ — 


+ S{N,z]mo)'- 


1 - - 

2 


(“A) 


1 - f (n + i) 


u 


mo 


du = 


V 


|m-mo| / 


27ri Jk{0,i) “ f (“ + ^) \/l - 2 :^ 

from Eqs. flll3p and fllldp we obtain, after simple calcnlations, Eqs. 


(114) 


(115) 
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